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Abstract 

The alpha and cluster decay properties of the ^^^"^^^Nd, '^^"^^^Gd, '^^"'''^Hg and 
Pb even-even isotopes in the two mass regions A = 130-158 and A = 180-198 are 
analysed using the Coulomb and Proximity Potential Model. On examining the clusters at 
corresponding points in the cold valleys (points with same A2) of the various isotopes of a 
particular nucleus we find that at certain mass numbers of the parent nuclei, the clusters 
emitted are getting shifted to the next lower atomic number. It is interesting to see that the 
change in clusters appears at those isotopes where a change in shape is occurring 
correspondingly. Such a change of clusters with shape change is studied for the first time in 
cluster decay. The alpha decay half lives of these nuclei are computed and these are 
compared with the available experimental alpha decay data. It is seen that the two are in good 
agreement. On making a comparison of the alpha half lives of the normal deformed and super 
deformed nuclei, it can be seen that the normal deformed Nd, " Hg and Pb nuclei are 
found to be better alpha emitters than the super deformed (in excited state) ' TS[d, ^ " Hg 
and Pb nuclei. The cluster decay studies reveal that as the atomic number of the parent 
nuclei increases the Nt^Z cluster emissions become equally or more probable than the N=Z 
emissions. On the whole the alpha and cluster emissions are more probable from the parents 
in the heavier mass region (A= 180- 198) than from the parents in the lighter mass region 
(A= 130-158). The effect of quadrupole (P2) and hexadecapole (P4) deformations of parent 
and fragments on half life times are also studied. 
Email address: drkpsanthosh@gmail.com 

1 Introduction 

The a-decay studies coupled with the cluster decay results have been used to point out 
the closed shell effects of the nuclei involved in the decay for quite some time now. The 
Z = 64 sub sheU was first predicted by Rasmussen et al. [1] from the study of a - decay 
energies and Schmidt-Ott et al. [2] noted it from the measurement of the a - decay reduced 



widths. Gupta and co-workers [3-6] have used the alpha and cluster decay results to predict 
shell effects in the daughter nuclei. Poenaru et al [7] predicted the spherical '"^"^Sn daughter 
radioactivity from their decay studies. Further a deformed radioactivity at Z = 74-76 and 
N = 98-104 was also noted by Gupta et al. [8] through cluster decay results. 

In this paper we have investigated the alpha and cluster decay properties of some 
even-even isotopes in the two mass regions A = 130-158 and A = 180-198. The nuclei 
which are studied in the present work are the ^^^"^^^Nd, ^^^"^^^Gd, ^^^"^'^^Hg and ^''^-wsp^ 
isotopes. All these isotopes are of particular interest for study because of the varied behaviour 
they present with respect to their shapes. The 150 < A < 190 mass region to which these 
isotopes belong is a well known region of deformation having various degrees of 
deformations and super deformations. There are a number of studies which throw light on the 
shape variations of all these neutron deficient isotopes. An abrupt change in the mean square 
charge radius, 5< r > was seen in Hg isotopes at mass numbers, A < 1 86 from isotope shift 
measurements [9-11]. The change in the mean square charge radius, 5< r > between A = 186 
and 185 was interpreted as a shape transition [12-14] while the isomeric shift in Hg was 
explained in terms of shape coexistence phenomena [15-17], a picture already suggested from 
nuclear spectroscopic measurements [18-21]. The even-even Hg isotopes with 180 < A < 188 
exhibit shape coexistence at low spins [19-33]. The coexistence of spherical and deformed 
degrees of freedom in light lead isotopes was proposed by a number of authors [17, 34-37]. 
The A = 144-158 mass region of the Gadolinium isotopes is a region in which deformation 
sets in rapidly and also the effects of the minor closed shell at Z = 64 for N < 88 have been 
observed. The even- even Gadolinium isotopes with 86 < N < 106 are shown to have positive 
ground state deformations in the Skyrme HF model [38]. Li the case of Gd the admixture 
of the negative deformation also seems to be important. 

The two mass regions under study also have the peculiarity that some of the nuclei 
among them are super deformed. Note that superdeformation here refers to the observation of 
super deformed (excited) bands in these nuclei, though their ground-state deformations are 
not very much different from other neighbouring nuclei. On the other hand, a deformed or 
normal deformed nucleus is the one where super deformed band(s) are not observed. Super 
deformed (excited) bands have been experimentally observed in " Nd, " Gd, " Hg 
and '^^'^^"^pb isotopes. In the present paper the alpha and cluster decay results have been 



analysed to get an idea of how the super deformed (in excited state) nuclei behave with 
respect to alpha and cluster decays as compared to their normal deformed counterparts. 

The Coulomb and Proximity Potential Model [39] have been used to study alpha and 
cluster radioactivity in various mass regions of the nuclear chart. Li this model the interacting 
barrier for the post scission region is taken as the sum of Coulomb and Proximity potential 
and for the overlap region a simple power law interpolation is used. The model (CPPM) was 
used to study cluster radioactivity and has predicted half life times for various proton rich 
parents with Z = 56-64 and N = 56-72, decaying to doubly magic ^ Sn. The model (CPPM) 
was also used to study the cold valleys in the radioactive decay of ^48-254^^ isotopes [40] and 
the computed alpha decay half lives are in close agreement with the experimental data. This 
model is modified [41] by incorporating the ground-state deformations P2 and P4 of the parent 
and daughter treating cluster as a sphere and the effect of deformation on the half-lives are 
studied, hi the present paper we have calculated the half lives and other characteristics of the 
chosen nuclei within the Coulomb and Proximity Potential Model. The effect of quadrupole 
and hexadecapole deformations of parents and fragments on the half lives is also studied. The 
details of the model are given in Sec. 2. The results obtained and the discussions therein are 
given in Sec. 3. The study of the effect of deformation on cluster decay half life time is given 
in Sec. 4. Finally the conclusions from our study are given in Sec. 5. 

2 The Coulomb and Proximity Potential Model 

In Coulomb and Proximity Potential Model (CPPM), potential energy barrier is taken 
as the sum of Coulomb potential, proximity potential and centrifugal potential for the 
touching configuration and for the separated fragments. For the pre-scission (overlap) region, 
simple power law interpolation as done by Shi and Swiatecki [42] is used. The inclusion of 
proximity potential reduces the height of the potential barrier, which closely agrees with the 
experimental result. The proximity potential was first used by Shi and Swiatecki [42] in an 
empirical manner and has been quiet extensively used over a decade by Gupta et al [43] in 
the preformed cluster model (PCM) which is based on pocket formula of Blocki et al [44] 
given as: 

(D(^) = _ - (5-2.54)' -0.0852(^-2.54)' fors<1.2511 (1) 

0{£) = -3A37Qiq)\—^] for s> 1.2511 (2) 
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Here O is the universal proximity potential. In the present model, another 
formulation of proximity potential [45] is used as given by equations (6) and (7). In this 
model cluster formation probability is taken as unity for all clusters irrespective of their 
masses, so the present model differs from PCM by a factor Po, the cluster formation 
probability. But we have included the contribution of both internal (overlap region) and 
external barrier for penetrability calculations. We would like to mention that within fission 
model [46-49] penetrability through the internal barrier (overlap region) represent the cluster 
formation probability. In the present model assault frequency, v is calculated for each parent- 
cluster combination which is associated with zero point vibration energy. But Shi and 
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Swiatecki [50] get v empirically, unrealistic values 10 for even A parent and 10 for odd A 
parent. 

The interacting potential barrier for a parent nucleus exhibiting exotic decay is given 



by 



V = ^^ + V^(z)+ ^ ^ ,forz>0 (3) 
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Here Zi and Z2 are the atomic numbers of the daughter and emitted cluster, 'z' is the 
distance between the near surfaces of the fragments, 'r' is the distance between fragment 
centres, i the angular momentum, // the reduced mass and Vp is the proximity potential 
given by Blocki et al [44] as 



V(z) = 47:yb 
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With the nuclear surface tension coefficient, 

r^0.95l7[l-lJS26(N-Zy/A^] MeV/fm^ (5) 

Here N, Z and A represent neutron, proton and mass number of parent, O represent 
the universal proximity potential [45] given as 

0(5) = -4.41e-"'°'"'' , for s>l. 9475 (6) 

0(£) = -1.7817 + 0.9270^7 + 0.0169 ^r'- 0.05148 £' for < s < 1.9475 (7) 

With s = z/b, where the width (diffuseness) of the nuclear surface b ~1 and Siissmann 
central radii Ci of fragments related to sharp radii Ri as 



C, = i?. 
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For Ri we use semi empirical formula in terms of mass number Ai as [44] 



i?, =1.284'" -0.76 + 0.84^'" (9) 

Using one dimensional WKB approximation, the barrier penetrability P is given as 



P = exp -^jV2/i(F-0Jz 



(10) 



Here the mass parameter is replaced by// = mA^A^ I A , where m is the nucleon mass 
and Ai, A2 are the mass numbers of daughter and emitted cluster respectively. The turning 
points "a" and "b" are determined from the equation V{a) = V{b) = Q . The above integral 
can be evaluated numerically or analytically, and the half life time is given by 

^ln2^ rin2^ 
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represent the number of assaults on the barrier per second 
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and A, the decay constant. Ey, the empirical zero point vibration energy is given as [51] 
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£,= 2^0.056 + 0.039exp 
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The Coulomb interaction between the two deformed and oriented nuclei taken from 
[52] with higher multipole deformations included [53, 54] is given as 
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where /?(,, =1.284"^ -0.76 + 0. 84 "^ • Here a, is the angle between the radius vector and 
symmetry axis of the /' nuclei (see figure 1 of ref. [53]). Note that the quadrupole interaction 
term proportional to y^jiAi is neglected because of its short-range character. 

3 Results and discussion 

We first present the alpha decay results for the ^^^"'^^Nd, ^^^"^^^Gd, '^^"^'^^Hg and 
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Pb isotopes and then analyse the cluster decay calculations. The driving potential 
(difference between interaction potential and Q value of reaction) of a compound nucleus is 
calculated for all possible cluster-daughter combinations as a function of mass and charge 
asymmetries. The Q values are computed using the experimental binding energies of Audi 
et al [55]. The driving potential versus A2 (mass of one fragment) graphs are plotted for each 



set of nuclei and are displayed in Figures 1-4. From the plots it can be seen that the potential 
energy minima occur at He and all the other typical clusters ( ' Be, ' ^ C, ^ ' ^ O, ' Ne 
etc.) for almost all the nuclei. One interesting result that comes up from these plots is the 
change of clusters at certain mass numbers of the parent nuclei, i.e., when we examine the 
clusters at corresponding points in the cold valleys (points with same A2) for the various 
isotopes of a particular parent nucleus we find that at certain mass numbers of the parent 
nuclei, the clusters emitted are getting shifted to the next lower atomic number. Thus in the 
cold valley plots (Figure 1) of the Gd isotopes it can be seen that the cluster Mg in the cold 
valley plot of Gd changes to Ne in the cold valley of Gd. Similarly Ar in the cold 
valley of ^"^"^Gd changes to ^^S in the vaUey of ^'^''Gd, again '^'^Ca becomes "^"^Ar and this goes 
on up to ^ Ge at Gd which is changed to ^ Zn in the cold valley of ^ Gd. The cold valley 
plots for Nd isotopes (Figure 2) show a similar trend with changes in clusters emerging at 
^^^Nd and ^^^Nd. The cluster ^^Mg in the plot of ^^^Nd changes to ^^Ne at ^^^Nd, then ^^Si goes 
on to become Mg and this continues till Zn at ^ Nd which is getting changed to Ni in 
the cold valley of ^^^Nd. Similarly the clusters ^°Ne, ^^S, ^^Ti, ^"^Fe and ^^Zn in the cold valley 
of '^^d are changed to ^°0, ^^Si, ^^Ca, ^''Cr and ^^i in the valley of ^^^Nd. Moving on to the 
cold valley plots (Figure 3) of the Hg isotopes changes in clusters emerges at ^^''Hg and 
^''^Hg. The minima at ^"^e, ^^S, ^^Ti, ^^Fe, ''^Zn and ^^Ge in the valley of ^^^Hg are getting 
shifted down to ^°0, ^^Si, ^^Ca, ^"Cr, ^^i and ^^Zn in the plot of ^^^Hg. In the cold valley for 
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Hg it is found that cluster changes are emerging at Ne, Mg, Fe and Zn i.e, Ne in 
the vaUey of ''^^Hg changes to ^^O in the plot of ^^^Hg. Similarly ^''Mg, ^°Fe and ™Zn in the 
valley of ''^^Hg changes to ^^e, ^°Cr and '°Ni in the valley of ^^^Hg. In the case of the Pb 
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parents there are changes at Pb but it is not a continuous one as can be seen at the other 
cases discussed above. 

The change in clusters appears at those nuclei where a change in shape is occurring 
correspondingly. As mentioned in the earlier paragraph the first change in clusters is 
observed at Gd. According to the tables of MoUer et al. [56] the quadrapole deformation 
parameters P2 of ^"^^Gd and ''^''Gd are -0.053 and respectively, thereby indicating a change 
in shape from oblate to spherical at Gd. For the low spin states in the neutron deficient Gd 
isotopes the deformation decreases with increasing neutron number and a spherical shell 
structure develops when the shell closure at N = 82 is approached. More precisely in ^ Gd a 
ground band is observed which has a structure resembling that of a vibrator or a rotor with a 
modest deformation. In contrast, Gd [57] show a structure which requires an explanation in 



terms of shell model excitations. This clearly indicates a shape change at the Gd isotope 
from deformed to spherical. Thus cluster changes are taking place at Gd where 
correspondingly a shape change is observed. Our studies show cluster changes at ^^"^Nd where 
there is a shape variation from normal deformed to superdeformed (in excited state). The 
Nd nucleus with experimental quadrupole deformation, P2 = 0.349 [58] is normal deformed 
prolate while the ^ Nd nucleus is superdeformed (in excited state). The next change in 
clusters is found at Nd. A shape change from prolate (P2 = 0.154 [56]) to oblate 
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(P2 = -0.138 [56]) is taking place when we go from Nd to Nd. Experimentally it is also 
found that Nd is superdeformed (in excited state) therefore there is also a change from 
superdeformed to normal deformed when going from Nd to Nd. Moving on to the Hg 
isotopes changes in clusters are observed at ^ Hg and Hg. In the case of ^ Hg we find 
shape variations from oblate to prolate i.e, ^Hg has P2 = -0.139 [56] while ^ Hg has 
p2expt.= 0.132 [58]. Similarly the oblate shape at ^''^Hg (p2= -0.130 [56]) changes to prolate at 
'^^Hg (P2expt. = 0.1 155 [58]). Also at '^^Hg there is a transition from superdeformed (in excited 
state) to normal deformed because ^^"^Hg has superdeformed excited bands while ^^^'^^^Hg 
nuclei are normal deformed. In the case of Pb isotopes no appreciable change in clusters is 
observed in our studies which mean there is no shape variation when going from Pb to 
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Pb. Recent experiments on the ground-state charge mean-square radii also show that the 
ground-state in all even-even Pb nuclei remains of spherical shape. Thus according to these 
observations, whenever a systematic change in clusters occurs at a particular isotope, a 
corresponding shape transition is occurring there. It is evident that our ground state cluster 
decay model is able to predict the shape change of nuclei in the ground state as well as in the 
excited state. 

3,1 Alpha decay results 

The logio(Ti/2) values corresponding to alpha decay of all the chosen Nd-Pb parent 
nuclei are plotted against the mass number of the parent and is presented in Figure 5. All the 
known experimental alpha decay half lives are included in the figure and it can be seen that 
they match well with our predictions. The experimental alpha decay data are from the 
refs. [59-68]. In the case of the Nd isotopes ^ Nd is the only nucleus around the chosen range 
which has experimental alpha decay data therefore its half life has been included in the graph. 
Along the same line the alpha decay calculations of ^^^Gd, ^^^Gd and I77.i79.i8i,i83,i85pjg ^^^^ 
been performed and included in the logio(Ti/2) graph for alpha decay. As far as the Pb 



isotopes are concerned the alpha half lives of all the nuclei (both even-even and even-odd) in 
the range A = 182-191 and again the alpha half lives of the even-even isotopes from A = 200 
to A = 210 have been estimated and presented in the alpha decay half lives graph. This has 
been done so because all the nuclei in the range A = 182-191 have experimental data while 
calculations for the isotopes from A = 200 to A = 210 have been conducted to include the 
lone experimental data of Pb. 

From the alpha decay half lives graph it is seen that ^^^ "^^^^Nd and i48.i49,150.i51.152q^ 
have their half lives within the present experimental limit of 10 s. Again with the exception 
of Hg all the other Hg nuclei considered (both even-even and even-odd) have their alpha 
decay half lives within the present experimental limit of 10 s. Similarly all the lead isotopes 
in the range 182 < Z < 200 (both even-even and even-odd) and the Pb nucleus have alpha 
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decay half lives within the order of 10 s. This means that all these nuclei are alpha instable 
and most of them have been experimentally detected as can be seen from the figure. Again 
from a first glance at figure 5 it is obvious that the elements in the heavier mass region 
(A = 180-198) are more prone to alpha decay than the ones in the lighter mass region 
(A =130-158). 

Figure 5 for the alpha decay results also presents a number of shell structure effects of 
both the parent and the daughter nuclei. The low T1/2 value for the ^'^'^Nd nucleus is due to the 
N = 82 closure in the Ce daughter. The extremely large half life of the Gd parent is 
attributed to the N = 82 spherical shell closure while the instability of 148,i49,150,15i,i52q^ 
parents are due to N = 82 and N = 82 shell closure in the corresponding Sm daughter nuclei. 
The lowest alpha decay half lives amongst all the nuclei are got for the ^ " Hg 
(Ti/2 = 0.059s for '^^Hg and T1/2 = 2.44xl0''s for ^^^Hg) and the i^2-i90p^ parent nuclei 
(Ti/2= 0.119s for ^^^Pb and Ti/2= 6.8xl0\ for ^>b). The daughter nuclei for the alpha decay 
of the ^'''"^^^Hg parents are the Pt isotopes with Z ~ 76 and N = 94, 96, 98, 100, 102, 104 and 
106 which implies the magic or semi -magic nature of these neutron shells. The N = 98, 100, 
102, 104 and 106 shell closures are again reflected in the Z ~ 82 Hg daughter nuclei for the 
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alpha decay of the " Pb parents. Several authors [3, 8, 69], have predicted Z = 76, 74 and 
N = 96, 98, 100, 102 and 104 as major (deformed) closed sheUs. Also the Z = 76 shell gap is 
got as a clear dip in the experimental shell energy curves [70] for Sr and Zr isotopes. Then 
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there is the very obvious extremely large rise in half life at the doubly magic Pb parent 
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nucleus (Z = 82 and N = 126) and again a sudden dip at the Pb parent which is of course 
due to the Z =^ 82 and N = 126 shell closures in the Hg daughter. 

On taking a comparison of the normal deformed and super deformed (in excited state) 
nuclei, it can be seen that except for the case of the Gd isotopes for all the other three set of 
nuclei, the normal deformed nuclei ( Nd, " Hg and Pb nuclei ) are better alpha 
emitters than the super deformed ones ( ' T^d, ^ " Hg and ^ Pb nuclei ). hi the case of 
the Gd isotopes the alpha decay half lives get influenced by the shell closure effects of the 
nuclei involved in the decay process. The super deformed ^"^^Gd nucleus is stable for alpha 
emission with a half life of 4.69x10 ^s while the super deformed ' '^^ Gd nuclei are alpha 
instable with half lives 2.05xl0^°s, 8.02xl0'^s and 9.09xl0^^s respectively. 

3.2 Cluster decay results 

The N=Z cluster and N^Z cluster results are studied separately for each set of parent 
nuclei. The logio(Ti/2) values for the N=Z clusters and Nt^Z clusters are plotted against the 
mass number of the parents and are displayed in Figures 6-13; the corresponding 
observations for each parent are presented below. The graphs for Nd parents show that there 
are almost equal number of N=Z and N^Z cluster emissions (six N=Z clusters and eight Nt^Z 
clusters) from them while in the case of the Gd parents there are more N^Z clusters being 
emitted than the N=Z clusters (five N=Z clusters and nine N^Z clusters). When we analyze 
the corresponding half life values we can see that for the Nd parents except for few emissions 
('^C, ^^O, ^^Si, ^°Si, ^^S from '^^Nd) almost all the cluster decays (both N=Z and N^Z) have 
their T1/2 > 10^ s and for the Gd parents except for two or three N=Z clusters ( C from 
^^°'^^^Gd and '^O from ^^^Gd) all the other cluster decays have their T1/2 > 10^°s. On 
comparing the cluster emissions from the Hg and Pb parents we find that there are equal 
numbers of N=Z and N^Z clusters (seven N=Z and seven Nt^Z clusters) getting emitted from 
Hg parents while there are more number of N^Z cluster emissions than N=Z emissions (three 
N=Z clusters and eleven Nt^Z clusters) from Pb parents. On making a study of the half lives 
of the cluster decays from Hg and Pb parents we can see that more than half the decays (both 
N=Z and N^Z) from the Hg and Pb parents have their T1/2 > 10^ s. Thus we may infer that as 
the N:Z ratio of the parent nuclei increases the number of Nt^Z clusters emitted are becoming 
equal to or greater than the number of N=Z clusters. Again from the analysis of the cluster 
decay half lives it is clear that the cluster decay rates for the parents in the heavier mass 
region are larger (smaller cluster decay half lives) than those for the parents in the lighter 



mass region. Thus the parents in the lighter mass region (A= 130-158) are found to be more 
stable against cluster decays than the parents in the heavier mass region (A=180-198). In the 
following sections we will discuss these results in detail for each set of parent nuclei 
separately. 

3.2.1 Nd parents 

Figures 6 and 7 give the logio(Ti/2) values plotted as a function of the parent mass 
number for the N = Z clusters and N ^ Z clusters respectively. The Nd parent is weakly 
stable against '^C, '^O (N = Z clusters) and ^'^Si, ^"^S (N ^ Z clusters) decays. Amongst these 
O cluster has the lowest half life of 1.68 x 10 s. The reason for the O decay from Nd 
being the most probable among the Nd parents is due to the proximity to the Z = 50 magic 
shell coupled with the midshell effect of the N = 64 neutron shell in the ^^ Te daughter. 
Interestingly the Si and S heavy cluster decays from Nd have their half lives 
(1.84 X 10 ^s and 2.96 x 10 s respectively) lower than the C light cluster decay 
(1.15 X 10^ s). It must be the stabilising effects of the N = 56 semi magic neutron number 
[70] which is the reason for the low half lives of the ^"^Si and ^"^S decays. The daughter nuclei 
for the ^°Si and ^"^S decays are the '°^Pd and ^^Ru nuclei with neutron numbers N = 56 and 
N = 54. 

The ^Be decay from ^^^"'^^Nd as well as the ^^C decay from ^^^"'^^Nd has their half 
lives, Ti/2 > 10 "^°s and therefore these cluster emissions must be considered as being hindered 
for the chosen Nd parents. Also from Figures 6 and 7 it can be seen that for the ^ TMd parent 
the ^ O, Si decay half lives are comparable with the alpha decay half life while for the 
^^^Nd parent the '^C, ^^O, ^^Ne, ^^Mg, ^Vg, ^°Si, ^^Si,^^Si and ^"S decays are more probable 
than the alpha decay. However the decay half lives in all these cases is very large, 
Ti/2 > 10^ s. Thus, although a number of shell effects are evident from the cluster decay 
results, the Nd parents are, on the whole, very stable against cluster decays. 

3.2.2 Gd parents 

The logio(Ti/2) results of Gd parents for N = Z clusters are displayed in Figure 8 and 
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those for N ^ Z clusters in Figure 9. The C cluster from Gd parent has the lowest half life 
(Ti/2=2.44 X 10 ^ s) among all the Gd parents. This low half life value of the C decay from 
Gd parent is due to the N = 82 spherical shell closure in the Ce daughter. Then there are 
dips for the Be and ^ O emission from the ^ Gd and ^^ Gd parents respectively. Be emission 
from ^ Gd isotope is due to the spherical N = 82 shell closure of the daughter Nd while the 



O decay from Gd is due to N = 82 spherical shell closure in the Ce daughter. The 
N = 82 spherical shell closure of the ^ Gd parent is obtained as a rise in T1/2 values for all 
clusters (except ^°Ne) at the '"^^Gd isotope. Of all these clusters only ^Be from ^^"^Gd and ^^C 
cluster from ^ Gd and ^^ Gd can be classified as weakly stable as all the other decays have 
relatively large half lives. 

Coming to the N ^ Z clusters, all the decays are either stable or very stable with the 
lowest half lives for the Si decay from Gd and the C decay from ^ Gd. The low value 
for the Si decay is due to the Z = 50 spherical shell closure in the Sn daughter while the 
C decay is bringing forth the N = 82 shell closure in the ^ Ce daughter which is also 
evident as a dip in the logio(Ti/2) versus mass number of parent graph. Another dip is got at 
the ^ Gd isotope for the O decay. This is again due to the N = 82 shell closure in the Ba 
daughter. The Sn daughter radioactivity is again coming up as a small dip at the ^^ Gd parent 
for the Si decay. It has also one of the lowest half lives (Ti/2=3. 33x10 s) among the N 7^ Z 
cluster decays. Although these decays have half lives much larger than the present 
experimental limit they are many orders smaller than the before mentioned N = Z cluster 
decays except for two or three cases. 

As can be seen from the logio(Ti/2) graphs there is a dip at the ^^"^Gd nucleus for a 
large number of clusters with the most probable one being the ^ O decay with a half life of 
8.04 X 10^ s. This decrease in half life values at the ^^ Gd nucleus is particularly visible in the 
case of the N 7^ Z clusters wherein it can be seen that the isotopes on the left and right side of 
the ^^ Gd nucleus are having larger T1/2 values than ^^ Gd for almost all the cluster emissions. 
Thus the '^"^Gd nucleus is found to be less stable against cluster decays as compared to its 
neighbouring isotopes. This decrease in half lives at ^ Gd is a classic case of daughter shell 
effects coming into play - almost all the dips in half lives at ^^ Gd is due to the shell closures 
in the corresponding daughter nuclei and they have already been discussed in the earlier 
paragraphs. 
3.2,3 Hg parents 

Figure 10 gives the logio(Ti/2) results for the N = Z clusters while Figure 1 1 gives the 

same for the N 7^ Z clusters. Of the probable N = Z clusters C from " Hg, Be from 
i76-i82pjg^ i6q g.^j^ i^^-isopjg^ 28g. g.^^^ i76pjg ^^^ 24^^ g.^^^ i76pjg parents havc half lives 
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within the present measurable limit of 10 s. Also the decay of C from Hg, Be from 
^^^Hg, "'O from '^^Hg and ^"^e from '^''Hg have half lives very close to 10^°s. Since all the so 



far observed daughters in decays of radioactive nuclei were magic or nearly magic nuclei the 
^^C emission from ^'''"^^^Hg, ^Be from ^^^"^^^Hg and ^^O from ^^''"^^^Hg indicate shell closures 
at '64-172^^ i68-i74Qg ^^^ i60-i64pj^ daughter nuclei. That is the Z = 76, 74, 72 proton shell 
closures and the N = 88, 90, 92, 94, 96, 98 neutron shell closures are coming up at these 
nuclei. This means predicting major closed shells at Z = 72 and N = 88, 90, 92 and 94 in 
addition to the already predicted Z = 76, 74 and N = 96-104 shell closures. Coming to the 
N ^ Z clusters, from figure 1 1 it can be seen that Si emission from Hg and Hg and 

0/^ 1 -7/' -5/^ 

Mg from Hg have half live values within the present experimental limit. The clusters Si 
and S from Hg also decay with half lives very close to the present experimental limit. 

As expected, the logio(Ti/2) versus mass number of parent graphs is showing dips at 
those masses whose decay involves daughters with known spherical magicities. Thus in the 
case of Hg parents, the minima occur at Hg parent for the decay of Si cluster with 
daughter ^ Dy having N=82 shell closure; then there are dips at ^^ Hg parent for Si cluster 
decay, ^^'^Hg parent for ^"^S cluster decay and '^^Hg parent for ^''S cluster decay, all of which 
have daughter nuclei ('"^^Dy, ''^*'Gd and ''^*'Gd respectively) withN=82 spherical shell closure. 
This is once again proof for the role of shell closure in cluster decay. 

The Hg and ^ Hg isotopes with just two protons less than the Z = 82 shell closure 
together with the proximity to the N = 102 and N = 108 neutron shell closures (both nuclei 
just two neutrons less than the N = 102 and N = 108 neutron shell closures respectively) are 
to be considered as stable nuclei. A strongly prolate neutron stable shell has been predicted at 
N = 102 by Hannachi et al, [33] while the shell closure at N = 108 has been emphasised by 
Bengtsson et al [70]. The analysis of the data on y-ray spectra by Carpenter et al, [71] has 
also revealed a prolate minimum for the ground state at N = 102. The decay results for the Hg 
isotopes present some interesting cases of instabilities at these stable nuclei - these are the 
^°Si and ^"^S decays from '^°Hg and the ^Be and '^C decays from '^''Hg. All these decays have 
half lives within or very close to the present experimental limit. The stability of the ^^ Dy 
daughter nucleus (with N ~ 82 and Z value exactly in the middle of two spherical magic 
numbers 50 and 82) is brought forth in the Si decay from ^ Hg while the S decay from 
Hg is due to the stability of the ^ Gd daughter nucleus. On a similar note the Z = 74, 76 

1 -J A 1 -TO 

and N = 100, 102 shell closures are again getting established from the W and Os 
deformed daughter radio activities in the Be and C decay from Hg. Gupta et al. [3] have 
observed similar cases of instabilities from ^ Ba and Hg nuclei. 



3,2,4 Pb parents 

Figure 12 which is the plot of logio(Ti/2) versus A, mass number of parent for N = Z 
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clusters shows that C cluster decay has the lowest half life among all cluster emissions from 
Pb parents with C from Pb having a half life of 5x10 s which is close to the present 
experimental limit of 10 s. The Be emission from ^ '^ Pb and ^ O emission from "^ Pb 
have their half lives in the range, 10 ^s < T1/2 < 10^ s and therefore they can be classified as 
being weakly stable. Coming to the N ^ Z clusters (Figure 13) Si cluster from Pb is 
weakly stable with the lowest half hfe (T1/2 = lO^^s). Also ^^Si from ^^^'^^b, ^^Ne and ^''S 
from '''¥b and ^^Mg from ^'^^'^^^Vh with half lives in the range, lO^^s < T1/2 < 10^°s are all 
weakly stable. 

The lowest T1/2 values of C emissions from ' Pb is attributed to the Z = 16 
proton shell closure coupled with the N = 104, 106 neutron closures in the ' Os daughter 
nuclei. Similarly the low half life value for the Be emission from Pb is due to the N = 106 
neutron shell closure in the ^ Pt daughter nucleus while the ^ O emissions from the 
I92,i94,i96p^ are due to the Z = 74 proton shell closure together with the N = 102, 104, 106 
neutron shell closures in the ' ' W daughter nuclei. 

Amongst the N ^ Z clusters the Si and the Si decays from ^ ' Pb parents have the 
lowest half lives. The daughter nuclei for the Si decays are the ' Er nuclei and the 
daughters for the Si decays are the ' Er nuclei respectively. Thus it can be seen that 
these emissions are occurring due to the N = 92, 94, 96 neutron shell closures coupled with 
the mid-shell effect of the Z = 68 proton shell in the Er isotopes. The Ne emission from 
^''^b has ^^'^Hf as daughter and ^^S emission from ^'^^Pb has ^^^Dy as the daughter. ^^°Hf has 
Z = 72 and N = 98 while ^ Dy has Z = 66 and N = 90 shell closures in them. Finally the 
^^Mg decay from ''^2,i94p^ j^^s the ^^^'^^Vb nuclei as the daughters. The N = 96, 98 neutron 
shell closures are being manifested in the ^ ' Yb daughter nuclei. Thus as in the case of the 
Hg parents the already predicted Z = 76, 74 and N = 96 -104 shell closures are again coming 
up at the Pb parents and in addition new shell closures at Z = 72 and N = 90, 92, 94, 106 are 
also coming up in the daughter nuclei. 
4 Effect of deformation 

In the present work the nuclear proximity potential for oriented and deformed (with 
higher multipole deformation) nuclei are done following the prescription of Gupta and co- 



workers [53] with universal proximity potential given in eqns. 6 and 7. The coulomb potential 
for the two deformed and oriented nuclei is computed using eqn. 13. In fission and cluster 
decay the fragments are strongly polarized due to nuclear force and accordingly their 
symmetry axes are aligned. In the present calculations we consider only pole to pole 
configuration. The proper inclusion of higher multipole deformations along with generalized 
orientation contributions may prove important in deciding the cluster decay paths of various 
clusters. 

Figures 14-18 and Table 1 represent the comparison of computed logarithm of half 
life times for the cluster emissions from '^^"^^^Nd, ^^^"^^^Gd, ^'^"^''^Hg and ^^^-^'^^'Ph parents for 
the cases of without deformation (spherical), with quadrupole deformation (P2) and with 
quadrupole and hexadecapole deformations (P2 & ^4). The experimental quadrupole 
deformations are taken from [58] and the rest from the tables of MoUer et al [56]. It is 
obvious from the table and figures that the half lives decrease with the inclusion of 
quadrupole deformation P2 due to the fact is that it reduces the height and width of the barrier 
(increases the barrier penetrability). We would like to mention that the sign of hexadecapole 
deformation have no infiuence on half life time. 
5 Conclusions 

An analysis of the alpha and cluster decay of the ^^^"^^^Gd, ^^^"^^^Nd, ^^^"^''"Hg and 

1 qo 1 go 

Pb even - even isotopes are carried out using the Coulomb and Proximity potential 
model. This set of parents comprise of both normal deformed and superdeformed (in excited 
state) nuclei. Therefore one aim of our study was to see how the superdeformed (in excited 
state) nuclei behave with respect to alpha and cluster decays as compared to their normal 
deformed counterparts. 

One important aspect of these nuclei is the shape variations they present as we move 
from one isotope to another and therefore their decay results were analyzed to see whether 
the shape changes infiuence their decay pattern in any manner. It was found that at certain 
mass numbers of the parent nuclei, the clusters emitted are getting changed from the ones that 
are emitted from the isotope just before it. More precisely, when we examine the clusters at 
corresponding points in the cold valleys (points with same A2) of the various isotopes of a 
particular parent nucleus we find that at certain mass numbers of the parent nuclei, the 
clusters emitted are getting shifted to the next lower atomic number and interestingly the 
change in clusters appears at some of those nuclei where a change in shape is occurring 
correspondingly. As already mentioned the first change in clusters is observed at Gd and 



correspondingly there is a change in shape from oblate to spherical at Gd. In the case of Nd 
isotopes the cluster changes are at Nd and ^ Nd. At ^ Nd there is a shape variation from 
normal deformed to superdeformed (in excited state) and a shape change from prolate to 
oblate is taking place at Nd. Here also we would like to mention a change from 
superdeformed (in excited state) to normal deformed at Nd. Finally in the case of Hg 
parents there are changes in clusters at ^ Hg and ^ Hg. At Hg we find shape variations 
from oblate to prolate and similarly at Hg the shape changes from oblate to prolate. Again 
we would like to mention that at ^ Hg there is a transition from superdeformed (in excited 
state) to normal deformed. In the case of Pb isotopes no appreciable change in clusters is 
observed in our studies which mean there is no shape variation when going from Pb to 
i>b. 

The effects of quadrupole and hexadecapole deformations of both parent and 
fragments on half life times are studied using Coulomb and proximity potential for oriented 
and deformed nuclei as the interacting barrier. It is found that inclusion of quadrupole 
deformation P2 reduces the height and shape of the barrier (increases barrier penetrability) 
and hence the half life time decreases. 

The most important part of our analysis was regarding the coming up of new shell 
closures (spherical or deformed) in either of the two mass regions (A = 130 - 158 and 
A = 180 - 198) under study. It is established that a large decay half life indicates that a parent 
nucleus is stable as far as alpha and cluster decay is concerned while a small half life is the 
indication for the daughter to be stable. Based on this idea we have come across a number of 
shell closures (some known and some new) in the different nuclei involved in the study. 

The elements in the heavier mass region (A = 180 - 198) are more prone to alpha 
decay than the ones in the lighter mass region (A= 130 - 158). The rise in the alpha decay 
half lives at the ^ Gd and Pb parents is attributed to the spherical sub-magic Z = 64, magic 
N = 82 and doubly magic Z = 82, N = 126 shell closures respectively. The N=82 spherical 
shell closure in Ce and Sm daughter nuclei are got as dips in the alpha half lives of the 
^^^Nd and ^^^Gd parent nuclei respectively. The normal deformed ^^^Nd, ^^^"^^^Hg and '^¥b 
nuclei are found to be better alpha emitters than the super deformed (in excited state) 

134, 136^^^ 190-194jjg ^^^ 194p^ ^^^j^- 

The cluster decay results show that on the whole cluster emissions are more probable 
from the parents in the heavier mass region (A= 180- 198) than from the parents in the lighter 
mass region (A= 130-158). The N = 82 spherical shell closure of the Gd parent is obtained 
as a rise in T1/2 values for all cluster emissions. The N = 82 spherical shell closure in the 



daughter nuclei is again getting manifested as dips in the cluster half lives for the cluster 
emissions from the Gd and Hg parents. The Be emission from ' Hg, C emission from 

^ "^ Hg, and O emission from ^' " Hg have their half lives within the present measurable 
limit of 10 s thereby indicating the possibility of major (deformed) closed shells at Z = 76, 
74, 72 and N = 88, 90, 92, 94, 96 and 98 in their ^^^-^^^W, ^"^"^^^Os and ^^^-^^^m daughter 
nuclei. Similarly the ^Be emission from ^^^''>b, the ^^C emission from i''2,i94p^ ^nd ^^O 
emission from ^ "^ ^Pb have the lowest half lives among all the cluster emissions from the Pb 
parents thereby bringing forth the Z = 76, 74 and N = 102, 104, 106 shell closures again. The 
Si and S decays from ' Hg isotopes present some interesting cases of instabilities in 
stable nuclei. The Hg and Hg isotopes with just two protons less than the Z = 82 shell 
closure together with the proximity to the N = 102 and N = 108 neutron shell closures are 
stable nuclei but the above mentioned decays from these nuclei have half lives within or very 
close to the present experimental limit. 
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Table 1. The comparison of computed values of logarithm of half life times for various 
clusters from ^^ "' Hg parents with and without including deformations. 



Parent 
nuclei 


Emitted 
cluster 


Daughter 
nuclei 


Q value 
(MeV) 


logio(Ti/2) 


Without 
deformation 


With 


With 
P2&P4 


_ '^"Hg 


%e 


l/2pt 


6.925 


-1.23 


-2.30 


-2.16 




«Be 


^^^Os 


13.298 


19.78 


17.29 


17.52 




16q 


l60Hf 


38.927 


21.99 


6.90 


9.45 




18q 


b8jjf 


31.312 


41.98 


27.37 


29.64 




^^Mg 


'^^Er 


62.683 


29.12 


6.31 


8.73 




^^i 


^^^Dy 


75.383 


29.23 


12.70 


12.74 




^^S 


^^^Gd 


69.161 


64.88 


68.12 


69.15 




^^S 


^^^Gd 


85.112 


33.70 


30.17 


34.34 




36g 


^^•^Gd 


80.444 


41.98 


15.32 


20.24 


'^^Hg 


%e 


lV4p^ 


6.578 


0.06 


-1.25 


-1.15 




^Be 


^^•^Os 


12.670 


22.16 


19.53 


20.02 




^^C 


'^% 


18.867 


45.15 


32.25 


31.83 




16q 


'^^Hf 


37.594 


24.63 


9.24 


12.55 




18q 


160Hf 


30.369 


44.63 


28.41 


30.39 




^^Mg 


'^•^Er 


56.666 


40.95 


19.76 


28.48 




^•^Si 


"*^Dy 
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Fig, 1, The driving potentials as a function of the mass number of light 
fragments (A2) for Gd isotopes. 
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Fig, 2, The same as for Fig.l but for Nd isotopes 
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Fig, 3. The same as for Fig.l but for Hg isotopes 
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Fig, 4, The same as for Fig. 1 but for Pb isotopes 
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Fig, 5. The logarithms of half lives for alpha emission from the chosen Nd, Gd, Hg and Pb 
parents calculated using CPPM and compared with the available experimental values. 
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Fig, 6. logio (T1/2) values for the N = Z clusters from Nd isotopes plotted as a 
frinction of the parent mass number. 
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Fig, 7. The same as for Fig.6 but for N ^ Z clusters from Nd isotopes. 
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Fig, 8. logio (T1/2) values for the N = Z clusters from Gd isotopes plotted as a 
function of the parent mass number. 
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Fig, 9. The same as for Fig.8 but for N ^ Z clusters from Gd isotopes. 
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Fig, 10. logio (T1/2) values for the N = Z clusters from Hg isotopes plotted as a 
function of the parent mass number. 
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Fig, 11, The same as for Fig. 10 but for N ^ Z clusters from Hg isotopes. 
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Fig, 12, logio (T1/2) values for the N = Z clusters from Pb isotopes plotted as a 
fiinction of the parent mass number. 
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Fig, 13. The same as for Fig. 10 but for N ^ Z clusters from Pb isotopes. 
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Fig, 14, Comparison of half life times for various clusters from '^^"'^^Nd parents with and 
without including deformations. 
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Fig, 15, The same as for Fig. 14 but for '^^"''"Gd parents 
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Fig. 16. The same as for Fig. 14 but for '"'"'"'Gd parents 
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Fig, 17, The same as for Fig. 14 but for '"^""''Hg parents 
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Fig, 18. The same as for Fig. 14 but for '"""''"Hg parents 
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Fig, 19, The same as for Fig. 14 but for '^^"'"'^Pb parents. 



